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Abstract 

We study a quantum-mechanical system of three particles in a one-dimensional 
box with two-particle harmonic interactions. The symmetry of the system is 
described by the point group Group theory greatly facilitates the ap¬ 
plication of perturbation theory and the Rayleigh-Ritz variational method. 
A great advantage is that every irreducible representation can be treated 
separately. Group theory enables us to predict the connection between the 
states for the small box length and large box length regimes of the system. 
We discuss the crossings and avoided crossings of the energy levels as well as 
other interesting features of the spectrum of the system. 
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1. Introduction 


During the last decades there has been great interest in the model of a 


harmonic oscillator conhned to boxes of different shapes and sizes 


2311 ■ Such 


model has been suitable for the study of several physical problems rangi^ 


from dynamical friction in star clusters 
and impurities in quantum dots 23|. 


^ to magnetic properties of solids [ 6 | 


One of the most widely studied models is given by a particle conhned to 
a box with impenetrable walls at —L /2 and L /2 bound by a linear force that 
produces a parabolic potential-energy function V{x) = k{x — xo)^/2, where 
|xo| < -h/2. When xq = 0 the problem is symmetric and the eigenfunctions 
are either even or odd; such symmetry is broken when xq 7 ^ 0. Although 
interesting in itself, this model is rather artihcial since the cause of the force 
is not specihed. It may, for example, arise from an inhnitely heavy particle 
clamped at Xq. For this reason we have recently studied the somewhat more 
interesting and realistic case in which the other particle also moves within the 


box 


24j |. Such a problem is conveniently discussed in terms of its symmetry 


point-group; for example: it is Ci when the two particles are different and 
C 2 h for identical ones. 

It follows from what was just said that the case of identical particles is 
of greater interest from the point of view of symmetry. For this reason in 
this paper we analyse the model of three interacting particles conhned to a 
one-dimensional box with impenetrable walls. In Section |2] we consider three 
particles in a general one-dimensional trap with two-particle interactions and 
discuss its symmetry as well as suitable coordinates for its treatment. In Sec¬ 
tion [3] we focus on the case that the trap is given by a box with impenetrable 
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walls and apply perturbation theory based on the exact results for inhnitely 
small box length. In Section 0] we discuss the limit of inhnite box length 
as well as the connection between both regimes. In Section |5] we construct 
a symmetry-adapted basis set for the application of the Rayleigh-Ritz vari¬ 
ational method and discuss the main features of the spectrum. Finally, in 
Section [6] we provide further comments on the main results of the paper and 
draw conclusions. 


2. Three particles in a one-dimensional trap 

We consider three structureless particles in a one-dimensional trap with 
a Hamiltonian of the form 

+ W{\Xi - X2\) + W(lX2 - Xsl) + W{\X3 - Xil), (1) 


where V(xt) conhnes each particle in a given space region (trap) and W{\xi — 
Xj\) are two-body interactions that couple the particles. It is convenient to 
dehne dimensionless coordinates {x,y,z) = (xi/L, X 2 /R, Xs/L), where L is a 
suitable length unit. The resulting dimensionless Hamiltonian is 


H' = 


2mL? 


H = - 


32 Q2 


dx‘^ 


+ 




dy"^ dz'^ 


+ v(x) + v(y) + v(z) 


+A [w{\x — y\) + w{\y — z\) + w{\z — x\ 


( 2 ) 


where n(a;) = 2mL‘^V{Lx)/h? axid Xw{\x — y\) = 2mL‘^W{L\xi — X 2 \)/h‘^, etc. 
In this equation A is a dimensionless parameter that measures the strength 
of the coupling interaction. 
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From now on we omit the prime and simply write H instead of H'. In 
addition to it, we assume that the trap is symmetric: v{—q) = v{q). The 
Hamiltonian Hq = H{X = 0) is invariant under the following transformations 



z) 


{x 


x,y, 

z) 


{- 

■x,y,z}p 


z) 


{x 

,-y,z]p 

x,y, 

z) 


{x 

,y,-z}p 

x,y, 

z) 


{- 

■X, -y,z}p 

x,y, 

z) 

-)■ 

{- 

■x,y,-z}p 

x,y, 

z) 

-)■ 

{x 

,-y,-z}p 

x,y, 

z) 


{- 

■X, -y, -z] 


( 3 ) 

where {a, &jc}p denotes the different permutations of the three real numbers 


a, b and c 


251. Note that the 


that is commonly named Oh 


8 coordinate transformations ([3]) form a group 
27|. In the particular case that v{q) = q^ 


m 


the symmetry of Hq is given by the full rotation group 0(3) 261. 

When A 7 ^ 0 the Hamiltonian is invariant under the 12 coordinate trans¬ 
formations 


{x,y,z) {x,y,z}p 

{x,y,z) {-x,-y,-z]p, (4) 


and we can describe the symmetry of the system in conhguration space by 


means of the point group D_ 


3d 


m 


27|. However, if v{q) = q^ the symmetry 


increases because the trapping potential is invariant under arbitrary rota¬ 
tions about any axis in conhguration space. The coordinate transformation 
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{x, y, z) = {qi, q 2 , qs) ■ J, where 


/ ^ L 0 

V2 V2 ^ 

1 1 _^ 

\/6 Ve \/6 

— — 1 

\ Vs Vs Vs j 

makes the coupling term independent of gs because 

, . f ^ V6 V2 V2 V6 " 

{x-y,y- z,z-x) = I V2gi, —ga - —gi, —^gi - —g 2 


( 5 ) 


( 6 ) 


Since the coupling term is therefore invariant under rotations about the gs 
axis by angles of j = 1, 2, 3,4, 5 the resulting Hamiltonian operator 


+A 


d2 Q2 q2 

+ 7 ^ + 


5gi2 5g22 5g32 

w (\/2 |gi|) + w 


I 2 I 2 I 2 

+ (?1 +^2 +(?3 


\/6 \/2 


+ w 


\/2 , \/6 

+ -Y^2 


( 7 ) 


is separable into its (gi, ga) and gs parts and exhibits symmetry D^h. 
If w(|s — f|) = {s — ty the resulting Hamiltonian 

/ Q2 q2 q2 


H = 


dqi"^ c^ga^ dq^'^ 
is fully separable and exhibits symmetry j27j|. 


+ gi^ + g 2 ^ + qs^ + 3A (ql + g^) 

n 


( 8 ) 


3. One-dimensional box 


As stated in the introduction we are interested in a system of three iden¬ 
tical particles in a one-dimensional box. To this end we resort to the Hamil¬ 
tonian ([2]) and choose the trapping potential 


n(g) 


0 if |g| < 1 
oo elsewhere 


( 9 ) 
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Therefore, the boundary conditions for the solutions to the Schrodinger equa¬ 
tion Hip = Eip become 




( 10 ) 


For concreteness we focus the discussion on the Hamiltonian 


H = - 


dx'^ 


d‘^ an 


[x-y) +{y-z) +{z-xy 


( 11 ) 


that is an extension of the two-particle model discussed earlier 


24| 


This problem is exactly solvable when A = 0 and the eigenvalues and 
eigenfunctions of Hq are 


TT 


^n\n2rLz 


{n\ + nl + Ug), ni, ^ 2 , ns = 1, 2,... 


^nin2n3 2/) 0ni(3^)0n2 (2/)0 713 ('^) 

A. i \ ■ + 1 ) 

pn[q) = sm---. 


( 12 ) 


All the sets of quantum numbers produced by the distinct permutations 
{ni,n 2 ,n 3 }p lead to the same energy (degenerate eigenstates). In addition 
to it, there is a Pythagorean degeneracy [2^ that we do not discuss here. 
Each of the functions cpn{q) is even (e) or odd (o) when n is either odd or 
even, respectively. A straightforward analysis based on group theory shows 


that 


25| 


{e,e,e}p 

-)■ 

Alg 

{e',e,e}p 

-)■ 

Mg 

ye'\e"'}p 


Mg 

{o,e,e}p 

-)■ 

Mu 




Em 
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{o,e'p 

e"}p 

-)■ 

Alu 

© 

^2u 

© 

Eu 

© 

Eu 

{e,o, 

, o}p 

-)■ 

Alg 

© 

E, 





(o, o', 

, e}p 

-)■ 

Alg 

© 

A2g 

© 

Eo 

© 

E, 

{ 0 , 0 , 

> o}p 

-)■ 

2^2u 







{o',o. 

> o}p 

-)■ 

2^2u 

© 

Eu 





(o, o', ( 

o"}p 

-)■ 

-^lu 

© 

7^2u 

© 

Eu 

© 

Eu, 


(13) 


where {s', s", s'"}p denotes a set of distinct permutations of products of func¬ 
tions 4>n{(l) with the indicated symmetry. The characters for every irrep, as 
well as some of the basis functions for them, are given in Table [1] 

It is instructive to carry out a straightforward calculation based on per¬ 
turbation theory of hrst order. For the hrst energy levels we have 

Svr^ 2 (tt^ — 6) 




EiEu 

E2Ai, 


^3Ai„ 


E2A2U 

E 3 A 2 U 


4 


TT^ 


-A + 0(A") 


1627r^ - 7297r2 - 4096 
37r2 1627r^ - 7297r2 + 2048 


A + 0(A2) 


2 

Ovr^ 


8l7r4 

2 (8l7r^ - 2437r2 


A + 0(A2) 


2048) 


8l7r4 


AE, - 


Ott^ 2 (8l7r^ - 2437r2 + 1024) 
~ ^ 81^^^ 
llTT^ 2 (Ovr^ - 38) 

4 Ott^ 

llTT^ 2 (Ovr^ - 38) 


A + 0(A2) 
A + 0(A2) 


^2Ea — 


Ott^ 


A + 0(A2) 
A + 0(A2) 


= 37r^ + 


- 3 


TT^ 


A + 0(A2) 


7e‘^ 1012507r^ - 275Q25tt^ - 2772992 


50 6 257r4 


A + 0(A2) 
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EzEu 

EiA\u 


Ttt^ 1012507r^ - 2 756 257r2 - 2048\/466441 


2 506257r4 

^ 1012507r^ - 2 756 257r2 + 2048\/46644l 
^ ^ 50 6 257r4 

77r2 , lOmOvr^ - 2 756 257r2 + 2772992 


+ 


506257r4 


A + 0(A2) 

A + 0(A2) 

A + 0(A2) (14) 


The eigenvalues that are degenerate at A = 0 exhibit the same parity u 
or g. Note that the energy levels -EsAig and E 2 Eg (stemming from the set 
of quantum numbers {3,l,l}p) remain degenerate at first order. We will 
discuss them in Section [5l 


4. Infinite box 

The parameter A introduced in Section [2] is proportional to the square of 
an arbitrary length L. In the particular case discussed in Section [3] the three 
particles are conhned into a box of length 2L as shown by the boundary 
conditions oa. Therefore, the inhnite-box limit L —)■ cxd corresponds to 
A —)■ cxo. We appreciate that when A —)■ cxd the motion of the system center 
of mass is unbounded and the eigenfunctions of the Hamiltonian operator 
are of the form 


/nin2fc(gi,g2,g3) = ^ni(gi)fe(g2)e*''®, (15) 

where gn{(l) is a harmonic-oscillator eigenfunction, ni,n 2 = 0 , 1 ,..., and 
—oo < k < oo. The eigenvalues behave asymptotically as 

lim X~^^'^Eni,n 2 ,k = 2\/3(ni + n 2 + 1). 

A—^CXD 


( 16 ) 










Since the symmetry of the eigenfunction is conserved as A increases we expect 
that the small-A eigenfunctions are connected with 


j cos{kq3)g ni (gi)fe(g2) 

[ sm{kq3)g 

ni (?l)fe(92) 


(17) 


instead of flT5|l . The symmetry of these functions is also given by equation 
m- For example, we expect that 


lim A 

A—^ oo 

lim 

A—^CXD 

lim X-^/^Eia 2 , 


lim A ^^‘^Eia 2 ^ = 2\/3 

A^cxd 

lim A-i/'Eie, = 4V^ 
lim X~^^‘^Eia^^ = SVS. 

X^oo 


(18) 


It is worth noting that {qi,q 2 ) and q^ are bases for the irreps E^ and A 2 U, 
respectively. 


5. Rayleigh-Ritz variational method with a symmetry-adapted ba¬ 
sis set 

As indicated in earlier pap ers on the application of group theory to a 


variety of problems 


24 


25 


29 


35| . we can obtain approximate eigenvalues 


and eigenfunctions of H from the eigenvalues and eigenvectors of the ma¬ 
trix representation of the Hamiltonian operator for every irrep S that 
we can treat separate from the other irreps. In the present case the linear 
combinations of eigenfunctions of Hq flT^ adapted to the symmetry species 


are 


Alg 


^2n-12n-12n-l 
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A 


2g 


(^2m-12n-12n-l + ‘f2n-l 2m-l 2n-l + V^2n-1 2n-l 2m-l) 

1 , 

{‘f2k-12m-12n-l + ^2n-l 2k-l 2m-l + ^2m-l 2n-l 2k-l 
+‘f2m-12k-12n-l + V^2n-1 2m-l 2fc-l + ^2k-l 2n-l 2m-l) 

(^2m—12n2n A ^2n2m—12n A ^2n2n2m—l) 

1 . 

W2k2m2n-l A ^2n-12k2m A ^2m2n-12k 
A(p2m2k2n-1 + A2n-12m2k + A2k2n-12m) : 


(19) 


{A2k-12m-12n-l + A2n-12k-l2m-l A '^2m-12n-12k-l 
—‘^2m-l2k-l2n-l — ‘^2n-l2m-l2k-l — 9^2fc-l 2n-l 2m-l) 


{'^2k2m2n-l A <^2n-12k2m A <^2m2n-12k 
A‘2m2k2n—1 A‘2n—12m2k A2k2n—12m) ) 


( 20 ) 


£ 


9 ■ 


(2v?2m-12n-12n-l ~ A2n-12m-12n-l ~ A2n-12n-12m-l) , 


1 

72 


{A2n-12m-12n-l ~ A2n-12n-12m-l) 


i2(p2k-l2m-12n-l — ‘^2n-12k-12m-l — ^2m-l2n-12k-l) , 


72 


{^2n-12k-12m-l — ^2m-12n-12k-l) 


i2cp2m-12k-12n-l — ^2n-12m-12k-l — ^2k-12n-12m-l) , 


10 



72 


(9^2n-12m-12fc-l ” 9^2fc-l 2n-l 2m-l) 


' 1 N 1 / 

i2(f2m-12n2n — ^2n2m-12n — ^2n2n2m-l) , (V^2n2m-12n ” ^2n2n2m-l) 


(2(p2k2m2n—l ^2n—\2k2m ^2m2n—12k') i iSP2n—12k2m V^2m2n—12fc) 


1 

.76 

' 1 N 1 

[2^2m2k2n-l ~ ^2n-12m2k ~ ^2k2n-12m) , W2n-12m2k ~ 9^2fc2n-12m)| , 

( 21 ) 


A 


lu 


1 . 

(V^2fc2m-12n-l + ^2n-12k2m-l + V^2m-1 2n-l 2A: 

~9^2m-12fc2n-l ~ </^2n-l 2m-l 2fc ~ V^2fc 2n-l 2m-l) 

1 


76 


{^2k2m2n + ^2n2k2m + ^2m2n2k 


^2m2k2n ^2n2m2k ^2k2n2m) i 


( 22 ) 


^2tj : 

^2n 2n 2n 

{^2m2n2n + ‘f2n2m2n + ^2n2n2m) 

1 . 

\SP2k2m2n H“ ^2n2k2m H“ ^2m2n2k 
-\-‘^2m2k2n + ‘^2n2m2k + ’^2k2n2m) 

(V^2m2n-12n-l + 9^2n-12m2n-l + 9^2n-l 2n-l 2m) 

1 . 

(V^2fc2m-12n-l + 9^2n-12fc2m-l + V^2m-12n-12fc 
+V^2m-12fc2n-l + 9^2n-12m-12fc + V^2fc2n-l 2m-l) , (23) 
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(2932m 2 n- 12 n-l ” 9 ^ 2 n- 12 m 2 n-l ” 9 ^ 2 n-l 2 n-l 2 m) , 


^Jc 2 (V^2n —1 2m 2n —1 9^2n—1 2n—1 2m) 


ye 


(2v32A:2m-12n-l ” 9^2n-12fc2m-l ” 9^2m-l 2n-l 2A:) , 


(9^2n-12fc2m-l ” 9^2m-l 2n-l 2fc) 


1 

Te 


(2932m-12fc2n-l ~ 9^2n-12m-12fc ~ 9^2fc2n-l 2m-l) , 


y2 


(9^2n-l 2ra—12k ^2k2n—12m—l) 


i^^2m2n2n ^2n2m2n ^2n2n2m) i 


V2 


iSP‘2n‘2m2n 9^2n2n2m) 


i^^2k2m2n ^2n2k2m ^2m2n2k') ) 


y2 


{‘^2n2k2m — ^2m2n2k] 


( 2932 m2fe2n — 9^2n2m2fc ” ^2k2n2m) , 


y2 


{‘^2n2m2k ~ ^2k2n2m) \ , 


(24) 


where k,m,n = 1,2,.... It is understood that <^ijk means that the three 
subscripts are different; equal subscripts are indicated explicitly as in (pm or 




*3 J- 


In this paper we have carried out a diagonalization of the matrix rep- 
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resentation of the Hamiltonian operator with 1000 basis functions of each 
symmetry species. Figs. [T]l6]show E{X) and for the hrst eigenval¬ 

ues of each irrep. The latter illustrate the conclusions drawn in Section 0] 
and in particular equations flT8|) . To facilitate the analysis we have drawn 
horizontal dashed lines that mark the limits fflB]) . 

It is well known that eigenvalues of the same symmetry do not cross but 
exhibit what is commonly known as avoided crossings 3^. Figs. [T]l6] exhibit 
some clear avoided crossings and others that appear to be actual crossings 
because of the scale of the hgures. As an example in Fig. [7] we show the third 
and fourth states of symmetry Aig on a hner scale to make it clear that they 
in fact undergo an avoided crossing. 

Fig. [H] shows a most interesting feature of the spectrum of this system of 
identical particles. The eigenvalues E 2 Eg and E^Aig are almost degenerate for 
small A (in Section |5] we have seen that they share the same energy of order 
zero and first order correction). However, the slope of the energy level E^^ig 
changes dramatically at the avoided crossing Ac with the level E^Aig but the 
slope of E 2 Eg does not change and the latter level approaches and remains 
close to for A > Ac. 


6. Conclusions 

It is interesting that a simple one-dimensional model exhibits such a rich 
symmetry. We have shown that group theory is remarkably useful for the 
interpretation of the results produced by approximate methods like pertur¬ 
bation theory and the variational method. The treatment of each irrep sep¬ 
arately from the others renders the analysis of the spectrum considerably 
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simpler. For example, the avoided crossings between levels of equal sym¬ 
metry can be studied separate from the crossings between levels of different 
symmetry. The interplay among the energy levels -£' 4 Aig and E 2 Eg is 

most interesting and one wonders if such behaviour also takes place in more 
realistic models of quantum-mechanical problems. 

Since the symmetry of the problem does not change with A we can resort 
to group theory to predict the connection of the states in the small-A and 
large-A regimes embodied in Eq. [TSl These results were conhrmed by the 
accurate Rayleigh-Ritz variational calculation in Section |5l 

Group theory is also useful for the calculation of matrix elements of ob¬ 
servables because we already know beforehand which of them vanish. This 
is the reason why we can treat each irrep separate from the others in the 
Rayleigh-Ritz calculation. The prediction of zero matrix elements is also 


the basis, 
transitions 


or example, of the well known selection rules in spectroscopic 


m, 


27|. 


In closing we mention that the combination of group theory and per¬ 
turbation theory has recently led us to the conclusion that the well known 
parity-time symmetry is less easily broken that other forms of the more gen¬ 
eral space-time symmetry 33(] . This result is of great relevance for the study 


of non-Hermitian Hamiltonians. 
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Table 1: Character table for D^d point group 



E 

2^3 

3^2 

i 

2^6 

ScTrf 


Mg 

1 

1 

1 

1 

1 

1 


Mg 

1 

1 

-1 

1 

1 

-1 


Eg 

2 

-1 

0 

2 

-1 

0 


Mu 

1 

1 

1 

-1 

-1 

-1 


Mu 

1 

1 

-1 

-1 

-1 

1 

X + y + z 

Eu 

2 

-1 

0 

-2 

1 

0 

{2z - X - y,x - y) 


+ y‘^ + 
xy + yz + zx 

{2z‘^ — x^ — t/^, x^ — t/^), 
{2yz — xy — xz, xy — xz) 



Figure 1: Eigenvalues for the irrep Aig 
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Figure 2: Eigenvalues for the irrep A2g 



Figure 3: Eigenvalues for the irrep Eg 
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Figure 4: Eigenvalues for the irrep Aiu 




Figure 5: Eigenvalues for the irrep A 2 U 
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Figure 6: Eigenvalues for the irrep Eu 
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Figure 7: Third and fourth eigenvalues of symmetry Aig 
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Figure 8: Third and fourth eigenvalues of symmetry Aig (blue, solid line) and second pair 
of eigenvalues of symmetry Eg (red, dashed line) 
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